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We report on a direct measurement of the mixing parameter y — (3.42 ± 1.39 ± 0.74)% in the 
D"— "D^ system by measuring the lifetime difference between the CP mixed final state K'^n~ 
and the CP even state K'^ K~ . We also present a study of the decay D'^ ^K^'k~ based on a 
sample of 149 ±31 observed events compared to 36 760 ±195 events observed in the Cabibbo 
favored channel D° ^K--k+. The observed branching ratio R = (0.404 ± 0.085 ± 0.025)% is 
used to obtain limits on the mixing parameters x' and y' and the doubly Cabibbo suppressed 
branching ratio, i?DCs. These studies are based on a large sample of photoproduced charm 
mesons from the FOCUS experiment at Fermilab (FNAL-E831). 



1 Introduction 



Mixing occurs because the neutral D mass eigenstates (or CP eigenstates in the limit of CP conser- 
vation) do not coincide with the flavor eigenstates and D . The mixing effects are parameterized 
by two dimensionless amplitudes 



AAf 



and y 



Ar 



(1) 



where AM is the mass difference between the two mass eigenstates, AF is the width difference and F 
is the average width. In the Standard Model, the D^ — D system mixing rateli?mix = + u"^)) is 
generally believed to be much smaller than the current experimental sensitivityl^. Nevertheless, recent 
measurements hint at a possible mixing effect at the edge of sensitivity. We report here on two such 
studies made with the FOCUS data. 

The data were collected by the FOCUS Collaboration during the 1996-97 Fermilab fixed target 
run in the Wideband Photon beam line using an upgraded version of the E687 spectrometer!. Charm 
particles are produced in the interaction of high energy photons {{E) ~ 180 GeV) with a segmented BeO 
target. In the target region, charged particles are tracked by 16 layers of silicon microstrip detectors 
which provide excellent vertex resolution. The momentum of the charged particles is determined by 
measuring their deflection in two oppositely polarized, large aperture dipole magnets with five stations 




Figure 1: (a) The decay reduced proper time distribution for the selected —*K~'k'^ and K'^ K~ events. The data is 
background subtracted and includes a small Monte Carlo correction, (b) The D*^ ^K^tv~ mass difference distribution 

with the signal and background fit contributions shown. 



of multiwire proportional chambers. Particle identification is determined by three multicell threshold 
Cerenkov detectors, electromagnetic calorimeters, and muon counters. 

2 Direct Measurement of AF from Lifetime Differences 

By measuring and comparing the lifetime for neutral D^s decaying to final states of pure even and odd 
CP a direct measurement of AF can be made. In this study!, the final state K^K^ is used as the 
CP even final state and, in the absence of a suitable CP odd candidate, the CP mixed state K~'k~^ 
is used. Assuming that K~'k^ is an even mixture of CP even and CP odd the relationship between 
the two lifetimes and the mixing parameter y is given by 

y _ even " ^CP odd _ t{D^K'k) ^ ^ 

even + TcPodd t{D^KK) ■ ^ ' 

The lifetime distributions for the 10 331 decays to K^K~ and 119 738 decays to i^^vr^ are shown 
in Fig. I^a. From exponential fits to the distributions we find t{D KK) = 395.7 ± 5.5 fs and 
t{D ^ Ktt) = 409.2 ± 1.3 fs where the errors are statistical (systematic errors are only evaluated on 
the ratio). Plugging these lifetimes into Eq. (^) we obtain y = (3.42 it 1.39 ± 0.74)%. 

3 Study of the Decay K+vr" 

The process IT may occur through either a doubly Cabibbo suppressed (DCS) decay, or by 

mixing to D followed by the Cabibbo favored (CF) decay D — >i^^7r~. The expected rate for the 
DCS decay relative to the CF decay (i?Dcs) is approximately tan'^^c — 0.25%, while Standard Model 
based predictions for i?mix range from 10~^ to 10~^. The large uncertainty in the D^ — D^ mixing 
rate is due to mixing mediated by intermediate meson states! whose contributions are not calculable 
in perturbative QCD. Nevertheless, large cancellations among the various intermediate meson states 
are expected and most studies conclude that the mixing rate should be at least a couple orders of 
magnitude less than 10~^ level S. Also, effects from beyond the Standard Model may enhance -Rmix- 

In this studyi, we begin by measuring the rate of — > K^-k~ decays relative to K~tt^ ^. 

The neutral D flavor is determined by requiring the decay chain D^* — s- 7:~^D^ — s- ir^^Kn). This is 

"Charge conjugate modes are implicitly included. 



achieved by looking for a narrow sig nal at 145 MeV/c^ in the mass difference between the D* and D 
candidates. The mass difference distribution for D^^K~^'k~ candidates is shown in Fig. ||b. We find 
149 lb 31 ^ K^7r~ events compared to 36 760 ± 195 decays of ^ K~7r~^ for a branching ratio of 
(0.404 ±0.085 ±0.025)%. 

The time dependence of the ^ K'^Tr~ decays is given by 



R{t) 



1 



(3) 



where x' and y' are phase rotations of x and y given by x' = x cos S + y sin 6 and y' = y cos 5 — x sin 6 
with 6 the strong force phase between the CF and DCS processes. Clearly, Eq. (P) indicates that in the 
case of significant mixing the measured branching ratio is dependent on the lifetime acceptance of the 
analysis. To account for this effect a large Monte Carlo sample of ^ K~7:~^ decays is used. Each 
Monte Carlo event accepted in the analysis is given a weight determined by the relative probability 
for an event with its lifetime given by Eq. divided by the probability for the same lifetime in the 
exponential decay rate used to generate the Monte Carlo. In this way we derive a relationship for 
-Rdcs as a function of x' and y' which depends only on the measured branching ratio, and the first and 
second moments of the accepted lifetime distribution in the Monte Carlo. The functional dependence 
on x' in the experimentally allowed region is small, while the dependence on y' (shown in Fig. ^ for 
the case of x' = 0) is large. For comparison, the y' and -Rdcs ranges from CLEO II. VQ and the y limit 
from FOCUS (discussed in Sect. |2|) are also shown in this figure. 

To determine limits on x' , y' and Rocs, the D — > Ktt data is split into high and low lifetime 
samples. The — > K^n^ branching ratio determined in each sample can be used to generate high 
and low lifetime curves in y' Rues space like the one shown in Fig. The point where these curves 
cross indicates the preferred values of y' and i?Dcs- To quantitatively determine the 95% confidence 
level (CL) allowed range, we integrate the likelihood for all points in the space and assign the 95% CL 
boundary to the high and low values beyond which the total integrated likelihood is equal to 2.5%. In 
determining y' (-Rdcs) limits the value of x' is set to zero (the value of x' with greatest likelihood) and 
we integrate over the entire allowed range .Rdcs (?/') variable. Using this procedure we find preliminary 
limits of 

-0.124 <y' < -0.006 and 
0.43% < i^Dcs < 1.73%. 

The large upper limit on .Rdcs and large negative lower limit on y' are the result of a second crossing 
point of the high and low lifetime curves. This property is an unfortunate side effect of using only one 
lifetime split - the fitting procedure! and limited statistics prevent more data splits. Nevertheless, 
the second crossing and its associated limits are far outside the allowed region of CLEO II. VQ and are 
also expected to be ruled out by the improved limits on x' and y' from FOCUS semileptonic mixing 
studies. 

To determine a limit on |x'| we integrate over the entire allowed range of y' and -Rdcs obtaining 
the preliminary limit of 

\x'\ < 0.039. 

We also determined a 95% CL boundary in x' y' space by integrating over all allowed values of i?Dcs and 
selecting the boundary line that is isometric in likelihood and contains 95% of the total likelihood. This 
boundary is shown in Fig. |2|b. Also shown in Fig are the best existing limits from the semileptonic 
mixing (E79i), lifetime differences (FOCUS) and D^^K+n- (CLEO II. V). 



4 Conclusions 



^0 



While the recent measurements in — D" mixing do not rise to the level of a discovery, they do 
warrant further study. A discovery of a non-zero y at the percent level (i?mix ~ 10~^) would not 




Figure 2: (a) The functional dependence of -Rdcs on y' with x' — 0. Fit values of y' and -Rdcs from CLEO II. Vu and 
the the FOCUS value for y from the lifetime difference study are shown for comparison. All errors are 1 a combined 
statistical and systematic, (b) The 95% confidence level allowed regions for x' and y' from the E791 semileptonic studyB, 
the CLEO II.V study of D°~^K+tv- and the FOCUS study of D°~*K+tv', and the FOCUS limits on y from lifetime 
differences. The direct measurement of y is only comparable to measurements on y' in the limit of strong phase 5 = 0. 
For the measurements of y' and y to be in agreement at the la level requires S > tv/AB. 



necessarily indicate new physics, but it would be unexpectedly large, and at the very least would lead 
to a deeper understanding of the processes involved in meson mediated mixing. 

The current data from FOCUS and CLEO II.V suggest two possible scenarios. The first scenario 
is that there is y-like mixing at the few percent level and a large (~ vr radians) strong phase. The 
second scenario is that some or all of the recent measurements are fluctuations. New measurements 
are required to determine which of these two possibilities is correct. 
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